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Abstract 

In this paper, a protocol called quantum energy distribution(QED) 
is proposed in which multi-parties can simultaneously extract positive 
energy on average from spin chains by use of common secret keys 
shared by an energy supplier. QED is robust against impersonation. 
An adversary, who does not have common secret keys and attempts 
to get energy, cannot obtain but give energy to spin chains. Total 
amount of energy transfer gives a lower bound of residual energy of a 
local cooling process by the energy supplier. 



1 Introduction 



Quantum teleportation [T] transfers any unknown quantum state to dis- 
tant places only by local operations and classical communication (LOCC). It 
has attracted much attention and been widely investigated. Today, it is con- 
sidered as a crucial building block of quantum communication. Recently, 
a new protocol named quantum energy teleportation (QET) in spin-chain 
systems has been proposed[3], which transports energy from one location to 
another only by LOCC. Entanglement of spin-chain ground states plays an 
essential role to realize QET. 

The protocol for QET [3J has been proposed for general spin chains with 
entangled ground states. Even before the advent of QET, spin chain systems 
have been hot topics of quantum information theory, because it is possible 
to apply it to short transmission of quantum information [4] . It is also known 
[S]that spin-chain entanglement is important to shed new light on compli- 
cated physical properties of ground states. 

In the QET protocol, a receiver of classical information from an energy 
supplier extracts positive energy from the ground state, accompanied by 
generation negative energy density in spin chain systems. Here, the zero of 
energy in the system is naturally defined by a value of the ground state. 
Though the concept of negative energy density is not so familiar to quantum 
information theory and quantum communication, it has been investigated 
in relativistic field theory for long time [6]. Detailed analysis for the spin 
chains can be seen in [3J. 

In this paper, an extended protocol is proposed, in which many authenti- 
cated consumers are able to simultaneously extract energy from the ground 
state by use of common secret keys shared by an energy supplier. Let us later 
call the protocol quantum energy distribution (QED). QED shows robust- 
ness against impersonation. Let us imagine that an illegal consumer appears, 
who does not have common secret keys and attempts to get energy from spin 
chains. Then we can conclude that the adversary does not obtain but give 
energy to the spin chains. We also notice that total amount of energy transfer 
of QED is related with local cooling. Local cooling is a short-time process in 
which energy is extracted from an excited system only by local operations at 
a certain site, without use of global time evolution generated by the system 
dynamics. In general, local cooling is unable to extract all energy of the ex- 
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cited system and residual energy remains in the system. The total amount of 
energy distributed via QED gives a lower bound of that residual energy of a 
supplier's local cooling for an excited state. We also analyze QET and QED 
protocols in the Ising spin chain system. Amount of energy transmission is 
evaluated depending on distance from the supplier. 

We confine our attention to short-time-scale processes in which dynamical 
evolution induced by the Hamiltonian is negligible. Meanwhile let us assume 
that classical communication between qubits can be repeated many times 
even in the short time interval. 

This paper is organized as follows. In section 2, we briefly review QET. 
In section 3, extending QET, a QED protocol is proposed. In section 4, we 
discuss a relation between QED and local cooling. In section 5, we analyze 
the Ising spin chain system and demonstrate the QED protocol. In the final 
section, conclusion is given. 



2 Brief Review of QET 



In this section, we shortly review QET. Detailed explanation is seen in 
[3J, including negative-energy physics of spin chains. Let us consider a very 
long spin chain system with Hamiltonian given by 

n 

where T n is the nth site energy density operator. In order to capture the 
essence of QET, let us focus on the nearest neighborhood interaction case. 
The operator T n is Hermitian and take the form of 

71+1 

t «=e n °™ ,7) ' 

7 m=n— 1 

where O^' 7 ' is a local Hermitian operator at site m. The ground state|g) is 
an eigenstate with the lowest eigenvalue of H. When we do not take account 
of gravitational interaction, absolute values of energy is irrelevant and just 
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difference of values makes sense. Hence, subtracting constants from energy 
density and the Hamiltonian, we obtain the following relations without loss 
of generality 

(g\T n \g) = 0, (1) 
H\g) = 0. (2) 
Due to Eq.(j2J), the Hamiltonian becomes nonnegative: 

H>0. 

In many models, \g) is a complicated entangled state. Using the entan- 
glement, Alice who stays at site ua can transport energy to Bob at site % 
only by LOCC. Taking account of the nearest neighborhood interactions, let 
us define localized energy operators of Alice and Bob as follows. 

n A +l 

Ha = /2 

n B +l 
n=nB — 1 

For later convenience, let us introduce several operators as follows. U A and 
U b are unitary Hermitian operators given by 

U A = n A -a nA , (3) 

U B =n B ■ a nB , 

where a are Pauli vector matrices, ua and n B are three-dimensional real unit 
vectors. The operator U A can be spectral decomposed into 

,u=0,l 

where Pa (/x) is a projective operator onto the eigenspace with an eigenvalue 
(— \y of Ua- Ub is time-derivative operator of Ub defined by 

U B = % [H B , U b ] = % [H, U B ] . 
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Next let us introduce two real coefficients as follows. 



t=(g\U B HU B \g)>0, 



(4) 




(5) 



cos (26) 



sin(2#) 



T) 



Finally we define a unitary matrix Vb ((J>) for \i = 0, 1 as follows. 



V B (/i) = / cos 9 + i (-If U B sin 0. 



(6) 



The parameter r\ is important for QET. If \g) is separable, we can generally 
prove that r\ vanishes. As seen below, QET transports no energy when r] = 0. 
Thus, in later discussion, we assume that \g) is an entangled state such that 



In order to perform QET, let us assume that Alice is a good distance 
from Bob such that 



and that Alice and Bob share many copies of spin chain systems in the 
ground state \g). Now let me explain the protocol explicitly. The protocol is 
composed of three steps as follows. 

(1) Alice performs a local projective measurement of the observable U a 
for the ground state \g). Assume that she obtains the measurement result \i. 
She must input energy Ea on average to the spin chain system in order to 
achieve that local measurement. 

(2) Alice announces to Bob the result /j, by a classical channel. 

(3) Bob performs a local unitary operation Vb (//) to his qubit at site n#, 
depending on the value of //. Bob obtains energy output E B on average from 
the spin chain system in this process. 



n A - n B \ > 5, 
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It is noticed [3] that the average input energy Ea is evaluated as 



E A = J2(9\PA^)HP A (fi)\g)>0. 



(7) 



M=0,1 

A positive amount of energy Eb is released to Bob's devices for the operation 
Vb (fj) in step (3). Eb is given by 



It is stressed that dissipation effect of transfered energy in channels can be 
completely neglected for QET because we transmit only classical information 
through a classical channel. 

In the above analysis, it has been argued that Bob actually obtains en- 
ergy from the spin chain system. However, even after the last step of the 
protocol, there exists energy Ea, that Alice first deposited to the spin chain 
by herself. Then a natural question arises. Does Bob extract positive en- 
ergy without any cost ? This apprarent paradox can be resolved from the 
viewpoint of entanglement breaking by Alice. Detailed explanation is seen in 
[3]. It is concluded that, based on a pledge of Ea, Bob knowing classical in- 
formation [i has borrowed Eb in advance from the spin chains. When global 
cooling induced by both long-time evolution of the system and extraction of 
energy makes the state approaching the ground state, the residual energy 
and negative energy — E B around site n B are compensated. 

In the QET protocol, classical channels for Alice to inform measurement 
results are not assumed private and secure. Therefore, anybody can extract 
energy from spin chains by listening to the measurement results announced 
by Alice. In the next section, an extended protocol is proposed in which 
legitimate multi-users can extract energy but illegal users are unable to steal 
energy from spin chains at all. 



3 Quantum Energy Distribution 



In this section, a QED protocol is proposed, in which M consumers C m urn — 
1 ~ M) can simultaneously extract energy from spin chains by use of secret 
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classical information sent by an energy supplier S. The protocol is an ex- 
tension of QET assisted by quantum key distribution(QKD). Let us consider 
that S stays at n — 0. Assume that the spin chain is so long that we are able 
to treat the number of sites as infinite and that the entangled ground state 
has a very large (or divergent) correlation length. Let us assume the sites of 
S and C m are separated from each other such that 

\ncj > 5 > 
\n Cm -n Cml \ > 5. 
Here let us introduce Us and Vc m as follows. 

U s = n s -a =J2 (-1)^5 

V m (//) = I cos 9 + i (—if U m sin 9, 

where Ps (/i) is a projective operator onto the eigensubspace with an eigen- 
value (-If of U s , 

U m W"m ' ®nc m i 

and fis and real normal vectors. The localized energy operators for 

those consumers are given by 

Hc m = ^ Tn - 
Let us also define a time- derivative operator of U m as 

U m = t[H } U m ]=i[H Cm , U m \. 

Consider that supplier S and any consumer C m share common secret 
short keys k for their identification, by which they are able to perform secure 
QKD in order for S to send secret classical information to those consumers. 
Because any protocol for QKD including BB84[7] is effective, we do not 
specify QKD protocols. Also assume that all C m and S share a set of many 
spin chains in the ground state Now let me explain the protocol explicitly. 
The protocol is composed of the following six steps. 

(1) 5* performs a local projective measurement of observable Us for the 
ground state \g). Assume that S obtains the measurement result fi. S must 
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input energy Eg on average to the spin chain in order to achieve that local 
measurement. E$ is evaluated as 

E s = J2(9\Ps(»)HP s (pL)\g). 

0=0,1 



(2) S authenticates C m by use of common secret short keys k. 

(3) S and authenticated C m 's generate and share sufficiently- long pseudo- 
random secret keys K via a protocol for QKD. 

(4) S encodes the measurement results // by use of K and sends to 
authenticated C m 's. 

(5) C m decodes the measurement reults \i by use of K. 

(6) C m 's perform a local unitary operation V m (fj) in Eq.Q to their qubits, 
depending on the value of \i. Each C m obtains energy output E m on average 
from the spin chains in this process. E m is given by 



F I 

m 2 



£ +r? 2 -£ 

Sm 1 /m Sr 



(9) 



where 

^ m = (g\u s u m \g)- 

After step (6), the quantum state is given by 

pqed= e (n^^j^^^^^fn^o*))- an) 

This QED protocol is robust against impersonation attack. Let us imagine 
that an illegal consumer Derick appears at site no, who does not have k and 
attempts to get energy from spin chains. Then we can conclude that Derick 
does not obtain but give energy to the spin chains. The reason is following. 
Because Derick cannot get no information about /i, Derick makes randomly 
two local operations Vd (0) and Vd (1) given by 

Vd (aO = I cos 6 + % (— Vf n D ■ a nD sin 6. 
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Then, instead of Eq. lflUI) . the final state becomes 
Pd = \Y, V ° V) f II ^ M ) p s (A*) \9) (9\Ps (**) ( II V l (/*) ) Vj (//) . 



Evaluation of the average localized energy around Derick is straightforward 
and gives a positive value such that 



Tr [p D H D ] = lYl M 00 ^Vb 00 l*> > °- 
Here we have used 



m'=o,i 



P 5 (/i) , Vj (//) ff D Vb (/i' 



and 

(g\Vl (//) ff D Vb (//) l<?> = <^|vj (a*0 #Vb (a*') !<?>■ 

Because the value of Tr [p D H D ] is positive, Derick must input energy on 
average to the spin chains without gain. 

Finally we add a comment that it is possible to array an infinite number 
of consumers in the most dense distribution by putting consumers at n = 5m 
for nonzero integer m. The total amount of energy gain by the consumers is 
defined by 



C 77 



2^ 

m^O 



4 Local Cooling by Energy Supplier 



In this section, we discuss a relation between QED and a local cooling process 
by the energy supplier S of QED. In step (1) of the previous QED protocol, 
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S must deposit energy Ea to the spin chain. Let us imagine that S stops the 
protocol soon after step (1) and attempts to completely withdraw Ea by local 
operations. By a similar argument in [3] , it is shown that this attempt never 
succeeds. In step (1), S breaks entanglement between 5"s qubit and other 
qubits and the entanglement cannot be recovered only by local operations. 
(Of course, for a long time interval beyond the short time scale that we 
have considered, local cooling is naturally expected to make residual energy 
approaching zero by an assist of dynamical evolution induced by nonlocal 
Hamiltonians. The time evolution is able to recover the entanglement broken 
by S.) Hence, there exists nonvanishing residual energy E r of the local 
cooling. Though explicit values of E r can be obtained for a special class 
of spin chain systems, including the Ising spin chain analyzed in the next 
section, the evaluation of E r is not so easy for general spin chains. However, 
Ec in Eq. ffTTT) generally gives a lower bound of E r . 

The reason is following. Let us consider a general local operation of S, 
which is expressed by use of //-dependent measurement operators M s (a,fi) 
satisfying 

J2Ml(a,fi)M s (a,fi) = I. 

a 

Then the quantum state after that local cooling by S is given by 

Pc = Y< M s(«> V)P S (/i) \g) (g\Ps 0) Ml(a, »). (12) 

The residual energy E r is evaluated as 



Er = j min Tt[ Pc H s ], (13) 



where H s is the energy density of S given by 



i 



The key point is that the value of E r can be calculated from the quantum 
state of QED. If S performs the above local cooling after the end of the QED 
protocol, the quantum state is transformed from that in Eq. ffTU]) to 



P ( qL = ^M s (a,n) J] V m (//) P s (ai) \g)(g\P s 00 J] v ™ M M s(®^)- 



Here it is easily proven that 



Tr 



Pqed h s 



Tr [p c H s 



because Ms (a, /x) and V m (//) commute with each other. Thus E r is rewritten 

as 



min Tr 

{M s (a lAt )} 



PqED S 



(14) 



It is stressed that the following relation should hold because of nonnegativity 
of H. 



Tr 



Tr 



PqedHs] + ^ Tr [pg D /ic„ 



> 0. 



(15) 



C ,7 



n 



-E, 



c- 



(16) 



Moreover, it is shown that 

E Tr [pQedH, 

where Pgko ^ s ^ e quantum state in Eq. ffTUl) . From Eq. ffTBI) and Eq. ffTBT) . we 
obtain 



Tr 



PQED H S 



-E c >0. 



This gives a relation what we want by taking account of Eq. ([T 



E r > E, 



c- 



(17) 

Thus it has been proven that Eq< in Eq. fllip gives a lower bound of E r . 

There may be a question whether the bound in Eq. f|T7|) is achievable or 
not. However, this is very nontrivial. One of the neccesary conditions is 
to achieve the equality of Eq. f|T5|) even if negative energy density appears 
in some region. Though the answer is not known for spin chain systems, 
equality of a similar relation does not hold for a free field in two dimensional 
spacetime [8]. Hence, it might be impossible to attain the bound in Eq. (jT7|) . 
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5 Ising Chain Analysis 



In this section, we demonstrate QET and QED protocols in Ising models 
without spontaneous symmetry breaking. Detailed properties of the model 
can be seen in [9] -[11]. Let us write the Hamiltonian as 

oo oo 

H = -h <- J Yl «^ - E 9i 

n=—oo n=—oo 

where h > J > and E g is a constant which shifts the eigenvalue of the 
ground state \g) to zero: 

H\g) = 0, 

and a z n and are Pauli matrices at site n. 

Let us introduce a parameter A = J/h(< 1). When A = 1, the system 
becomes the critical Ising model. The energy density operator at site n is 
defined by 

T n = ~ha z n - -< « +1 + a x n _ x ) - e, 

where e is a real constant to satisfy (g\T n \g) = 0. The Hamiltonian can be 
expressed as a sum of T n : 

n 

The system can be mapped into a Fermionic system and solved analytically [9]. 
It is shown that correlation functions are evaluated, for exmaple, as 

(g\<\g) = (gK\g) = o, 

(gK\g) = G(0), 
(sKK' + ^i) \ 9 ) = 2G(-1), 

(gW v m ° y m+n\g) = A ( n ) 

G(l) G(0) ■■■ G(2-n) 
G{2) G(l) ■■■ G(3-n) 

G(n) G(n-1) ■■■ G(l) 
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where the function G{n) is defined by 



G(n) = L(n) + \L(n + 1), 

LM = i f - cos(fcn) dk. 
n Jo VI + A 2 + 2A cos k 

A relation 

(g\< K + ^-i) l^) = o 

is also obtained explicitly. For the case with A = 1, G(n) is calculated 

G(n) = 2 -t±L. 

Then the correlation functions are given by 

2 

(gWn\g) = -, 

7T 

(9\oZ(ot + o*_ 1 )\g) = ± 

7T 



(gWm<Tm+n\g) 



2 y 2 2n(n-i) h ( n y 



where 



7T 



n-l 



(4n 2 - 1) /i(2n) ' 



h(n) = Y[k n ~ k . 



k=l 

The asymptotic behavior of A(n) for large n is given by 

A(n ~ oo) ~ -le 1 ^ 1 / 12 ^ 3 ™- 9 / 4 , 

where the constant c is evaluated as c ~ 1.28 [TTj . 

Now let us first consider QET. To specify the protocol, we set 

Ua = o* a , 

and 

Ub = < b - 

The energy input Ea by the energy supplier of Eq. (J7J) is evaluated as 
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E A = -(g\oi A Hai A \g) = h(g\a z nA \g)+J(gK A « +1 + < A „ X ) \g) = hG(0)+2JG(-l) 



where we have used a relation as 



Ps^) = \[l + {-lf< A ]- 
For the case with A = 1, E A is given by 

E A = -h. 

7T 

The coefficient £ in Eq.(j4]) is evaluated as 

£ = {g\o* B Ho* nB \g) = 2h(g\a* nB \g) = 2«7(0). 
The time-derivative operator of Ub is given by 

U B = i[-ha" nB , al B ]=2hal B . 
Hence, the value of r\ in Eq.Q is calculated as 

V = ( 9 \U a Ub\9) = 2h(g\al A al B \g) = 2hA(\n A - n B \). 
From these values of £ and 77, the energy output E B of QET is computed as 



2h 

HiB — 



7T 



l+(^A(|n A -n B |)) 2 -l 



It should be stressed that nonvanishing values of E B are obtained for general 
values of A, including noncritical models. It is noted that the evaluation of 
E B is simple for the critical Ising case with A = 1. The asymptotic value of 
E B with A = 1 is obtained for \n B — n A \ ~ 00 from Eq. (l2Ui) as follows. 



E B ~-[-A(\n B -n A \) 

7T VZ 



7T 



-9/2 



(23) 
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Next let us consider QED with an infinite number of consumers in the 
most dense distribution. When A = 1, the total amount of energy transfer 
Ec in Eq. ffTTT) can be evaluated explicitly as 



E, 



2h 



c 



TV 



£ 

m^O 



I (|A(5|m;)V I 



6.2 x 10" 5 /i. 



(24) 



As discussed in section 4, this value gives a lower bound of E r of local cooling 
by S. In this solvable model, we can check explicitly the relation in Eq.( fT71) . 
The minimization of E r in Eq. ([T3|) among local operations is possible. 
For general values of A, the localized energy H$ is expicitly written as 



J {atx + of) ~ 3e. 



H s = ~ha z - Jo'l [a\ + a x _ x ) 

It is noted that the following relation holds for p c in Eq. ([12 
Tr 

= (g\ 



I ^x —X —X j I 'z i _z\ o , 



(-^K - ^-i°-2 ~ J + *?) - 3c) \g). 

By use of the above relation, we are able to manipulate as follows. 
Tr [p c H s ] = Tr[p c (-ha z -J^(af + a^))] 

+(9\ ("^2 " ^-1^-2 " J {oli + *1) - 3e) \g) 
= (g\ {hoi + J°o K + <ry ) \g) - Tr [p c {ha* + Ja* (a? + o^))] + (g\H s \g). 
By substituting (g\Hs\g) = and 

(<?| {hoi + K + ^-i)) \g) = hG(0) + 2JG(-1), 
it is obtained that 



Tr [ Pc H s ] = hG(0) + 2JG(-1) - Tr [p c (ha* + Jag (a? + cx^))] 
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Here it is useful to write Ps (//) as follows. 

P s (fi) = \fi,n = 0) (p, n = 0\l[l n = ~[I+ (-If <} 

After messy calculations using Eq. ( |25|) . we obtain 

Tr [ Pc H s ] = hG(0) + 2JG{-1) - hTr [p s a z ] . 
In general, the following inequality holds. 

1 > Tr [ Ps al] . 
The equality is attained by measurement operators as 



(25) 



(26) 



M A (/i = 0) 
M A (p = 1) 

without the a degree of freedom. These lead to the final result as followed. 



E r = h(G(0) - 1) + 2JG(-r 



When A = 1, the value of E r is given by 
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7T 



1 J ~ 0.91J. 



Because evaluation of the total amount of energy transfer Eq in Eq.l fTTl) is 
difficult, comparison between Ec and E r is not performed easily. However, 
by taking A = 1, the comparison is possible as follows. From Eq. ( f24j) . 
Eq. lfTTl) is verified because 0.91 > 6.2 x 1CT 5 . In this case, the bound of 
Eq. lfTTl) is not so tight. This result depends on the choice of Ua and Ub in 
Eq. lpTI) and Eq. (l22l) . If we change the choice, the bound may become tighter. 
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6 Conclusion and Discussions 



In this paper, a protocol for QED is proposed, in which many consumers 
can simultaneously extract energy from spin chains by use of common secret 
keys shared by an energy supplier. In this protocol, what consumers need for 
energy gain is just classical information about measurement results. Hence, 
dissipation process in energy transportation via channels can be neglected. 
This protocol for QED is robust against impersonation. An adversary, who 
does not have common secret keys and attempts to get energy, cannot obtain 
but give energy to the spin chains. We have also pointed out that the total 
amount of energy distributed via QED gives a lower bound of residual energy 
in a supplier's local cooling for a state excited by the supplier's measurement. 
Finally, QET and QED protocols have been studied in the Ising spin chain 
model. Amount of energy transmission is explicitly evaluated depending on 
distance from the supplier. Finally, the lower bound of Eq. fll7l) has been 
explicitly checked. 

For practical situations for QED, there remain some open problems of 
QED. They are listed below. 

One of them is related with energy dissipation. It is expected that dissi- 
pation effects in the energy transport of QED are severely suppressed even if 
a zero-temperature uncontrolled environment is coupled with the spin chain. 
What those consumers need for energy gain is just classical information about 
the measurement result without receiving energy directly from the supplier. 
This aspects is quite a contrast to ordinary energy transmission in the spin 
chains. In the transportation, excitations as energy carriers in the spin chain 
get gradually annihilated dependent on environment interaction properties. 
Detailed analysis of comparison between QED and ordinary energy trans- 
ports in spin chains is an interesting problem. However, the analysis is out 
of scope of this paper and will be discuss elsewhere. The finite temperature 
effect is also considered to be of importance for QED. 

In the QED protocol, we fix a unitary operation of B which takes a form 
in Eq. ([6]). However, the optimal local operation of B to extract maximum 
energy from the spin chain is not obtained yet. The optimal operation is 
crucial for analyses of local cooling because the extracted energy gives a 
lower bound of residual energy of local cooling by the supplier. 

Realistic implemetation proposals of QED are not reported yet. However, 
there are possible candidates. One of them might be the carbon nanotube. 
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Carbon nanotubes are able to contain fullerenes with atoms inside. The 
fullerenes acquire spins by doping some atoms. By laying the fullerenes side- 
by-side in a nanotube, a spin chain might be constructed and useful to check 
QED. Quantum dots and SQUID qubit systems might allow to create spin 
chains for QED. 
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